Does EELS haunt your photoemission measurements? 
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It has been argued in a recent paper by R. Joynt (R. Joynt, Science 284, 777 (1999)) that in the case 
of poorly conducting solids the photoemission spectrum close to the Fermi Energy may be strongly 
influenced by extrinsic loss processes similar to those occurring in High Resolution Electron Energy 
Loss Spectroscopy (HR-EELS), thereby obscuring information concerning the density of states or 
one electron Green's function sought for. In this paper we present a number of arguments, both 
theoretical and experimental, that demonstrate that energy loss processes occurring once the electron 
is outside the solid, contribute only weakly to the spectrum and can in most cases be either neglected 
or treated as a weak structureless background. 



I. INTRODUCTION 

Photoemission has been used for years as a reliable 
technique for probing the the electronic structure of the 
occupied states in solids ranging from insulators through 
semiconductors and metals to superconductors. In his 
article Joynt provided very convincing and interesting 
arguments that especially for badly conducting samples 
(roughly po > Q.lmVtcm, the Mott value) the photoelec- 
tron spectrum may be affected by energy loss structures 
resulting from the interaction with the time dependent 
fields set up by the photoelectron receding from the sur- 
face of the solid. He argued that the influence of these 
loss processes can be so strong that the spectrum is dom- 
inated by them and that therefore the intrinsic informa- 
tion regarding the electronic structure of the solid all but 
disappears. Since photoemission is playing such a promi- 
nent role in the discussion of strongly correlated materi- 
als like the High Tc's, or more generally the transition 
metal oxides, as well as in Kondo and Heavy Fermion 
systems, it is of quite some importance to further inves- 
tigate Joynt 's assertions. In this paper we study Joynt 's 
arguments and provide both experimental and theoret- 
ical findings that show that the effects due the losses 
discussed by Joynt are only a small contribution to the 
total spectrum and the zero energy loss probability for 
the photoelectrons dominates for samples of either good 
or bad conductivity. 



II. INTRINSIC AND PSEUDO-INTRINSIC 
EFFECTS 

Compared to other techniques, photoemission provides 
the most easy and direct measurement of the one elec- 
tron Green's function and the directly related occupied 
density of states of a solid, if one keeps in mind the 
possible influence on the photoelectron of a number of 



'pseudo-intrinsic' effects ||. By this we mean first of 
all the matrix element needed for the description of the 
amplitude of the optical transition probability from an 
occupied state to a high energy unoccupied state in the 
solid, followed by energy loss processes occurring on the 
electron's path to the surface and finally the description 
of the escape of the electron through the surface region 
into the vacuum. Where we emphasize that in the so 
called 'sudden approximation' we neglect all interactions 
and interference effects between the high energy excited 
electron and the hole left behind. 

Also one has to contemplate the truly intrinsic effects: 
in a one electron approximation the associated one hole 
Green's function is a delta peak at an energy determined 
by the band dispersion of the occupied states. In real- 
ity the electrons in the solid are usually not simple free 
electrons but they interact with other electrons, phonons, 
magnons etc. , resulting in one electron Green's functions 
now including a frequency and momentum dependent 
self-energy. For weakly interacting systems the initial 
delta function spectrum for such an electron broadens 
(asymmetrically) and attains a frequency distribution for 
each momentum vector, which basically provides infor- 
mation not only of the quasi-particle dispersion and life- 
time, but also of the way the electron is dressed inside 
the solid, due to the response of its environment to its 
presence or absence. In strongly interacting systems this 
self-energy causes a rather large spreading out of the ini- 
tial delta peaks describing the one hole Green's function, 
and the description in terms of a quasi-particle with a cer- 
tain lifetime may break down completely. In these cases 
it is indeed difficult to separate the intrinsic properties of 
the one hole Green's function from the pseudo-intrinsic 
effects due to the energy losses suffered by the excited 
electron on its way out of the solid. These losses are ba- 
sically dominated by the self-energy of the excited elec- 
tron. It is extremely important therefore to have good 
estimates of the contributions due to energy loss pro- 
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cesses to the spectrum so that these can be identified, 
and possibly corrected for, if they are substantial. 

Experimentally there are several ways of checking the ex- 
pected influence of these energy loss processes. The most 
obvious is to study the energy loss spectrum of electrons 
incident on the solid with an initial kinetic energy equal 
to that of the escaping photoelectron. These loss spec- 
tra provide one with direct information about the self- 
energy of an excited electron in the unoccupied states 
of the solid. For high energy electrons (Efc > 60 keV) 
a transmission EELS measurement is possible in which 
the obtained energy loss information is mostly due to 
bulk processes. However, to achieve high energy resolu- 
tion, photoemission measurements are usually performed 
at low photon energies (Eph < 100 eV) such as Hel radi- 
ation (21.22 eV) or even less. At these low energies elec- 
tron energy loss processes can be studied with reflection 
High Resolution EELS but one must then realize that 
the incident electron hardly penetrates the solid surface 
and that therefore the surface loss function is measured 
rather than the bulk losses experienced by a photoelec- 
tron originating from inside the solid as in photoemis- 
sion. The surface and bulk loss function are linked to 
each other however and thus information about the in- 
trinsic processes can be retrieved also from a reflection 
experiment, although the relative intensities of certain 
bulk losses compared to their surface analogs may differ. 
We note here that Joynt concentrates on losses occurring 
after the photoelectron has left the solid and so these 
should be directly related to the reflection ELS losses. 

Another way of getting the self-energy information of 
the excited electron is to study the spectral function in 
a photoemission measurement of a narrow atomic core 
level of the solid at a photon energy such that the ex- 
cited electron will have the required (low) kinetic energy. 
Although this loss spectrum will be intertwined with the 
satellite structure due to the self-energy effects involved 
in the sudden creation of the core hole itself, they will 
nonetheless provide us with an upper limit on the im- 
portance of the materials energy loss contributions to a 
photoemission spectrum. This was in fact suggested and 
used in Ref. j| to argue that the broad and intense en- 
ergy distribution seen in angular resolved photoelectron 
spectroscopy of the high Tc's was not a result of energy 
loss processes but mainly a direct result of the strongly 
energy dependent self-energies in these strongly corre- 
lated materials. A more detailed study of this effect in 
the High Tc's has recently been published ||. 

Thus, in the interpretation of spectra, it is generally as- 
sumed that the photoelectron intensity distribution is a 
true reflection of the single particle spectral function, and 
that the aforementioned pseudo-intrinsic losses can be 
identified and reckoned with if necessary. 



III. EXTRINSIC EFFECTS 

We now will focus on extrinsic broadening of photoe- 
mission structures such as the Fermi edge. Within pho- 
toemission this up to now meant just the finite instru- 
mental resolution, but using the picture that Joynt put 
forward in his recent article, we now also have to consider 
losses suffered by the photoelectron after it has left the 
solid and is on its way to the detector. These losses are 
directly comparable to the loss spectrum of a reflection 
EELS measurement. Joynt argues that these losses more 
than anything will severely distort any sharp feature, 
such as the Fermi edge, especially in the case of badly 
conducting solids and anisotropic materials. He claims, 
depending on the properties of the material under study, 
that this extrinsic distortion can be so dramatic that in- 
stead of observing a sharp Fermi cutoff in the spectrum, 
the observed spectral distribution will look like that ex- 
pected for a material with a so called pseudogap at Fip. 
If correct this would make photoemission unsuitable for 
the study of the one electron Green's function of badly 
conducting solids such as many of the colossal magne- 
toresistance materials and the High Tc's, among many 
others, are. 

Joynt substantiates his statement by deriving an expres- 
sion for the energy loss probability of an electron once 
the electron has emerged from the surface. He calcu- 
lates, in a classical picture, the average energy lost due 
to the interaction of the electron's time dependent elec- 
tric field acting on the polarizable metal left behind. The 
response of the metal is approximated by a Drudelike be- 
haviour. He then distributes this average energy lost over 
an energy loss spectrum using a probability distribution 
as a weighting function. The basic assumption is that the 
Born approximation is valid so that each electron suffers 
at most a single scattering event. This results in: 



2e 2 CL{u) 
nhvuj 2 



(1) 



in which P{ui) stands for the probability that the electron 
loses an energy u>, C is a constant for which Joynt gives 
a value of ~ 2.57, and L(u>) represents the loss function 
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We agree to a large extend with this derivation except 
perhaps for the constant C on which we will focus our 
attention later on in this paper. Let us for the moment 
assume that the Born approximation and Joynt 's deriva- 
tion is correct. We see that the material related prop- 
erties enter the equation for the energy loss distribution 
via the frequency dependent dielectric constant of the 
material. Which is not unexpected since this describes 
the response of the material to the time dependent field 
produced by the outgoing electron. Let us, as Joynt did, 
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take as an example the dielectric function e(u>) given by 
the Drude model: 



e{u>) = cr(uj) = — ; r 

y ' lu y ' w(l - iujr) 



(3) 



and plot P(ui) for different values of the resistivity po — 
1/cto and scattering time r (Fig. 1) one sees immediately 
that although the weight of P(w) is distributed differ- 
ently in each curve, the integral J °° P(ui)dcj reaches the 
same value in the end. There is indeed a well known 
general sum rule Q connected to this formula which is 
independent of material parameters: 
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This is a general sum rule that holds for any model of 
e(w) provided that it's a causal function. It depends only 
on the velocity of the outgoing photoelectron which is 
presumed constant in the process, consistent with the 
Born approximation valid for weak scattering. So, for 
an electron with a kinetic energy of 20 eV, as used in 
Joynt's examples, the integral from zero to infinity gives 
w 0.0843 for the total of the losses. 

There is of course another sum rule which conserves the 
integrated intensity of the spectrum. In this sum rule we 
must include the finite probability Po that an electron 
suffers no loss at all, so that 



l = Pn 



P(u)du 



(5) 



Using the results above we find that Po « 0.957, and 
therefore the normal single particle DOS will dominate 
the photoemission spectrum. 

Up to now we've just been concentrating on the free elec- 
tron part of the response function, of course there are 
more contributions (such as phonons, and interband tran- 
sitions) that are contained in the total dielectric function 
of a real material. Therefore the calculations presented 
by us and by Joynt, using only the Drude model, will al- 
ways overestimate the influence of the free electron con- 
tribution. 

If we take a simple example to illustrate this: our first 
sum rule pins down the total amount of losses from zero 
to infinite frequency, this means that other processes such 
as phonons, not contained in the Drude model will just 
steal away some of the weight carried by the free electron 
losses that we have considered so far. This implies that in 
a good metal, where the phonon part will be nearly fully 
screened by the free electron excitations, the (low en- 
ergy) loss spectrum is in first approximation indeed well 
described by the Drude model. If the material is then 
changed into a bad conductor the phonon part of the loss 
spectrum will gain more and more strength in the low en- 
ergy region of the spectrum at the expense of the Drude 



part because of the sum rule. In both cases however, Po 
will have a fixed value. Furthermore, any phonon contri- 
bution will not influence the PES spectrum in a smooth 
way on the low energy scale around the Fermi energy, but 
rather produce a step in the convoluted photoemission- 
EELS spectrum, since these phonon losses are peaked at 
the phonon frequency, and are never overdamped. The 
same holds for a good metal where there is a clearly de- 
fined loss peak at the plasma frequency. 
Thus, we can conclude that at 20 eV, J °° P(uj)dui « 
0.0843 is an upper limit for the losses due to only the 
Drude part. This very general result is in direct con- 
flict with the assumptions made by Joynt that Po could, 
within the approximations made by him, be very small, 
as he hints at the fact that processes other than the free 
electron losses will further reduce Po, and therefore he 
takes Po to be a fit parameter. In our opinion it is not 
possible to make an independent choice for the value of 
Po as Joynt did since it is in essence determined by the 
sum rules and is independent of material constants. 
From an experimental point of view, we know from re- 
flection EELS experiments (see e.g., Fig. 4) at incoming 
energies of around 20 eV, that the elastic peak (which 
represents Po) is by no means close to zero for any ma- 
terial. Only for very low incoming energies (below w 10 
eV) or when special surface waveguidelike conditions are 
met H can the zero loss peak be strongly suppressed. Be- 
sides this, in reflection EELS P(w) is twice as strong as 
the loss probability in this photoemission scenario, since 
the electron can lose energy both on the incomming and 
on the outgoing trajectory. In fact, one can use the same 
method as that used in calculating the reflection EELS 
loss probability (see Ibach & Mills ||) for this photoe- 
mission problem. It is interesting to note that one gets 
the same result except for a different numerical factor, 
[j. This presumably stems from a difference in Fourier 
transform convention and in Mills's case, avoiding in- 
tegrals such as equation (1) in Joynt's paper, which is 
not readily solvable analytically. Using the prefactor ob- 
tained by the procedure described by Mills the equation 
for the losses reads: 



P(u) = - — Im 
nujv 



-1 



1 + e(w) 



(G) 



Which then for the sum rule means that the losses are 
in fact much more severe than with Joynt's original pref- 
actor: we now have J °° P(ui)dcu w 0.65 again at a ki- 
netic energy of the electron of 20 eV, leaving Po at 0.35. 
Even here, the zero loss part will still be large enough to 
dominate the Fermi cutoff, if one thinks in terms of the 
complete dielectric function being involved. 
This result does imply however that working in the Born 
approximation is no longer valid, and a strong interac- 
tion picture containing also multiple losses needs to be 
applied, which is less straightforward to derive for a con- 
tinuous spectrum of excitations. In the case of discrete, 
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welldefined plasmon losses due to core level excitations, 
the loss spectrum has in fact been calculated by Langreth 
H with the result that multiple plasmon losses are seen 
distributing the energy loss over a wider energy range at 
the expense of the low energy losses. On the basis of this 
calculation we argue that multiple scattering corrections 
will not strongly influence the sum rules but instead they 
will just redistribute the losses over a wider energy range 
thereby reducing their influence in the low energy loss 
region. 

If we for the moment stick to the single scattering sce- 
nario and take for the losses only the Drude contribu- 
tion, we can calculate the photoemission spectrum for 
the same parameters as Joynt used, except now taking 
Po = 0.35. This is depicted in Fig. 2, the lower panel 
also contains the original calculation from Joynt, using 
Po = 0.01. From this we see, that although the lineshape 
of the photoemission spectrum is affected by the losses in 
the case of the bad metal (where the Drude contribution 
is overestimated!), there is still finite weight at the Fermi 
Energy and therefore the effect does not create a clear 
pseudogap structure. 

IV. EXPERIMENT 

As a last discussion point, we will present the case of 
Lai.2Sri. s Mn 2 07, a double layered colossal magnetore- 
sistance oxide with a ferromagnetic metal (at low tem- 
perature) to paramagnetic insulator transition at 125 K 
grown by the traveling solvent floating zone method ||. 
This material is a good candidate for testing Joynt 's as- 
sumption that materials with high resistivity will be more 
prone to the influence of losses on the Fermi region in 
photoemission spectra than those with low resistivity, as 
the resistivity changes by roughly two orders of magni- 
tude from below to just above the phase transition (see 
Fig.3) |§. We performed both reflection EELS at 20.5 
eV incoming electron energy and angle integrated PES 
using a Hel source. For both measurements the sample 
was cleaved in situ at a temperature of 60K, and at a base 
pressure of 8 x 10~ n in the case of EELS and 4 x 10~ n 
in the case of PES. As these samples deteriorate even 
at these pressures in a matter of hours, we ensured that 
measurements were performed within 2 hours after the 
cleave, before the peak at 9 eV binding energy in the 
PES spectra started appearing, which is associated with 
a change in oxygen stoichiometry at the surface jLl[ . 
The satellite and background corrected regions around 
Eip of the photoemission spectra are shown in Fig. 5 for 
T=60K (solid), 140K (dashed) and T=180K (dotted). 
The inset depicts the full spectra, taken at 60K before 
and after the temperature cycle. 

Since we performed EELS at a finite incoming angle with 
respect to the surface normal (6*™ = 35 deg) we have to 
apply a correction factor as described by Ibach and Mills 



(|] to extract the surface loss function Im (—1/1 + e(w)), 
and then multiply this by e 2 /hcuv to get P(lo) as de- 
scribed in equation (||) in order to use it to simulate a 
PES spectrum. Therefore, in Fig. 4, top panel are de- 
pict the EELS spectra as taken at 50K (solid) and 150K 
(dashed) including the zero loss line, and in the lower 
panel P(w) calculated from the data after subtraction of 
the zero loss line. 

We then use this P(ui) to calculate its influence on a 
PES spectrum assuming a constant density of states, and 
using various values of Po- This is depicted in the lower 
panel of Fig. 5. It can be seen from this figure, that 
unless Po = the picture of Joynt doesn't reproduce the 
photoemission spectrum at all. For any finite value of Po 
combined with a finite DOS at the Fermi energy, one will 
get a finite Fermi cutoff. To get a rough estimate of P 
from experiment, we can integrate the zero loss peak sep- 
arately and compare it to half the integral of the entire 
loss region up to 10 eV as measured in our EELS spec- 
tra, provided of course that we keep in mind that we used 
our detector in a mode which selects electrons within an 
opening angle of 2° around the specular reflection and is 
therefore not fully angle integrated, which makes us un- 
derestimate the losses relative to the zero loss probability. 
However, if we proceed in this way, we get for both the 50 
and 150K a ratio of P : P(w) = 0.82 : 0.18 which shows 
at least that Po is not close to zero, and neither does Po in 
our experiments depend on temperature (or resistivity) 
of the sample. Our findings agree with ARPES measure- 
ments by Dessau and Saitoh et al. |l^,|ll| in which they 
use Joynt 's argument that he expects the loss effect to 
be angle independent to show that therefore the angle 
at which the smallest change is observed in going from 
above to below Tc is indicative of the maximum magni- 
tude of the effect, and in their experiments turns out to 
be negligible. 

V. CONCLUSIONS 

In conclusion, we have argued that Joynt indeed raises 
an important question regarding the influences of extrin- 
sic losses on low energy photoelectrons, but we disagree 
with the statement that the losses will take the upper 
hand in determining the shape of the spectrum around 
the Fermi energy as there is a sum rule that renders the 
zero loss probability substantially finite. This at least 
holds down to photon energies such as the often used Hel 
line (21.22 eV), but may become a point of concern when 
really low photon energies are used. Of course, since the 
losses are by no means a small perturbation in this classi- 
cal, single scattering approach a full quantum mechanical 
treatment including multiple scattering is called for. Fur- 
thermore, we are not able to reproduce Joynt's formula 
exactly as far as the pre-factor is concerned, we believe 
however that the treatment by Mills (?]] is self-consistent 
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and avoids integrals with questionable convergence. We 
also have shown in the case of Lai^Sri S Mn 2 7 that 
we cannot reproduce the photoemission spectrum using 
a finite density of states up to the Fermi energy together 
with a finite value for P , so that we must conclude that 
there is a true pseudogap in this material both below and 
above the phase transition. 
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FIG. 1. Top panel: The loss probability P(u>) calculated for 
different values of resistivity p and relaxation time r, keeping 
the surface plasmon energy uj% p = 47r/(pr(l + too)) fixed at 
0.75 eV. Lower panel: Demonstration of the sumrule for the 
loss probability P(o>). 
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FIG. 2. Comparison between Joynt's original calculation 
with P = 0.01 and / P(w)duj = 0.0843 (lower panel) and 
our calculation using Mills' prefactor for the loss probability 
P(a>) and obeying the sum rule Po = 1 — / P(u>)du> = 0.35 
(top panel) The parameters are specified in the figure and T 
= 38K. Although the shape of the spectrum is affected, there 
is no clear sign of a pseudogap for the bad metal in our case. 
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FIG. 3. The resistivity of our Lai.2Sri.sMn20r samples 
in the ab plane (circles) and the c direction (squares) as a 
function of temperature (from Ref. 9). 
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FIG. 4. Top panel: Reflection EELS spectra, as taken, of 
Lai, 2 Sri,gMn 2 07, at 50K (solid) and 150K (dotted). Incom- 
ing/outgoing angle is 35° with respect to the surface normal. 
Incoming electron energy is 20.5 eV. Zero Loss FWHM = 60 
meV. The filled black curve is the 50K spectrum divided by 
150 to show the relative weakness of the loss features with 
respect to the zero loss electrons. Lower panel: Geometry 
corrected spectra, according to Ref. 6 multiplied by e 2 /huiv 
to get P(w). Inset: first 400 meV, showing the stronger pres- 
ence of a surface phonon in the insulating regime, relative to 
the metallic phase. 
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FIG. 5. Top panel: PES spectra at different temperatures 
above and below the phase transition in La1.2SV1.8M/i2O7 
close to the 'Fermi energy', the inset shows the spectra up to 
12 eV binding energy, before and after the temperature cy- 
cle ensuring the sample has not degraded in the mean time. 
Lower panel: Calculated PES spectra using the loss proba- 
bility P(w) constructed from the EELS data of Fig. 4 and 
assuming a constant DOS, for various values of the zero loss 
probability Po as indicated in the figure. 
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